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Using ideas from BCS descriptions of systems of fermions, a covariant extension of the relativistic 
Fermi gas model is presented as a way to incorporate correlation effects in nuclei. The model is 
developed for the BCS nuclear ground state and for final states consisting of a single plane-wave 
nucleon plus a BCS recoiling daughter nucleus. The nuclear spectral function is obtained and 
from it the superscaling function for use in treating high-energy quasielastic electroweak processes. 

QQ , Interestingly, this model has the capability to yield the asymmetric tail seen in the experimental 

f^ ' scaling function. 

o: 

p^ , PACS numbers: 25.30.-c, 25.30.Fj, 24.10.-i, 24.10.Jv 

I. INTRODUCTION 

in 

^^ ■ The issue of su per scaling in electroweak interactions with atomic nuclei has been the focus in recent studies [l|, y, 

y, 3, [a, [a, 0, d, [3, [l3, [ill Ej Ej [iJ, m, M, Il3| because of its intrinsic interest and because it offers a way to obtain 

reliable charge-changing and neutral-current neutrino-nucleus cross sections. The latter are usually essential when 

addressing fundamental neutrino properties, specifically neutrino masses and the neutrino oscillations that result from 

those masses. A summary of the basic concepts involved in superscaling analyses is provided below in Section HIl 

^ ' Owing to the complexity of nuclear dynamics it is not obvious that the nuclear response to an electroweak field 

(~| \ superscales; indeed for this to occur the single-nucleon cross section should factorize, which is almost never completely 

'""'■ true. Certainly in extreme on-shell independent-particle models such as the Relativistic Fermi Gas (RFG) this is valid, 

-^1 \ although in general factorization is seen to be violated to some degree, since the nucleons in the nucleus are not on- 



O 






►^ , shell and the independent-particle model is only an approximation. For instance, in other studies done within the 
pg ' context of the shell model (for instance, [l^, [ll|) one sees breaking of factorization from off-shell effects, typically 
(^\ , at the few percent level for the kinematics of interest when applying the ideas of superscaling, namely, momentum 
CN ■ transfers of order 1 GeV/c and energy transfers near the quasielastic peak position. Clearly, if collective effects (for 
'"^ \ instance, giants resonances) are important, the independent-particle description breaks down and thus superscaling 
f — ■ should not be expected to hold in the low momentum and energy transfer domain where such effects occur. Finally, 
^^ ' when meson-exchange currents play a significant role the superscaling violations that may or may not occur require 
00 ; detailed study [ll [la US UlI ■ 

^^ ■ In the present work we shall not address these general issues, but instead shall focus on an extension of the RFG 

^ , approach, namely, we shall include only on-shell nucleons in an independent-particle model where factorization is exact 
and put our efforts into going beyond the degenerate description provided by the extreme RFG model. To do this we 
resort to a model inspired by the BCS theory of condensed matter physics (see, e.g., [22|) in which electron Cooper 
J3 ' pairs are correlated via an attractive interaction mediated by the lattice. In nuclei correlations among pairs of nucleons 
also occur 23] , in this case mediated by NN interactions including both the long-range attractive residual interaction 
(the pairing force) and also the strong repulsive force that occurs at short distances. However, and remarkably, the 
BCS wave function can be shown to hold valid (in the thermodynamic limit) in both instances in the sense of being 
an eigenstate of the well-known pairing Hamiltonian independently of the sign of the interaction (of course in the 
attractive case the BCS state provides a minimum for the system's energy, in the repulsive one it yields instead a 
maximum) . 

Hence we shall deal with the problem of assessing the impact of correlations on the superscaling function employing 
the BCS formalism with appropriate modifications, such as retention of covariance, to adapt it to the high-energy 
physics of atomic nuclei. This is addressed in Section Hill after having briefly outlined in Section |TT] the basics of the 
superscaling approach for the longitudinal nuclear response i?i, namely the part of the total inclusive electroweak 
response that is believed to superscale the best. 

The numerical results obtained with our approach are presented in Section HVl To summarize briefly our findings, 
we shall see that, likely because in the BCS spirit we limit ourselves to an independent quasi-particle description of 
nuclear matter, scaling of the first kind (independence of momentum transfers q) appears to occur not only at the 
quasielastic peak (QEP), but also at both lower and higher energy transfers lu. We find that the onset of first-kind 



scaling already occurs at momentum transfers of order 500 MeV/c, whereas away from the QEP the onset only 
occurs at quite large momentum transfers (of the order of 2 GeV/c). Furthermore, interestingly the shape of the 
superscaling function turns out to be non-symmetric around the QEP, being larger to the right and smaller to the left 
of it, namely, in agreement with experiment and thus lending support to our approach. However, in our model when 
in the so-called scaling region (below the QEP) first-kind scaling is reached as a function of q from below, which is not 
what is experimentally found. Finally, scaling of the second kind (independence of nuclear species) is also explored 
in Section [TVl and shown to be relatively well satisfied, given that an appropriate momentum scale is chosen for each 
nuclear species. 

In the Conclusions (Section |V| we briefly discuss the significance and limitations of our model. 

II. THE LONGITUDINAL RESPONSE FUNCTION 

As anticipated in the Introduction we shall concentrate on the longitudinal EM channel; all other elec- 
troweak responses can be developed using similar arguments to those presented in the following. We start 
by briefly summarizingthe basics of the supersca ling approach, which has been discussed in depth in previous 
work [l|, 0, y, |j, [^, [y, 0, [a, M, UQ, ^M, Il3- llSl. Uj, [H, \l^- Here we employ the same notation and formalism as 
in the paper by Cenni-Donnelly-Molinari [3\ (to be referred to as CDM). 

Within the framework of the Plane- Wave Impulse Approximation (PWIA) the longitudinal response function of a 
nucleus to an external electroweak field bringing three-momentum q and energy lo into the system reads 

RL{q,cj) - // dpd£p'^2nS{p,£)^^nL{q,u;,p,£), (1) 

JJs PQ Ep En 

where Ep = ^p^ -|- m^ and E^ = \/(p + q)^ + 'm'j^ are the energies of the initial (struck) and final (ejected) 
nucleons, respectively, with m^ the nucleon mass. Indeed in PWIA it is assumed that only one vector boson is 
exchanged between the probe and the nucleus and that this one is absorbed by a single nucleon. The probability of 
finding one nucleon in the system is provided by the spectral function of the system, S{p,£), which depends on the 
missing momentum p and an energy £, the latter being essentially the missing energy minus the separation energy and 
given explicitly by the excitation energy of the residual nucleus in the reference frame where it moves with momentum 
—p. The factor TZl in Eq. ([1]) is the longitudinal single-nucleon response which is in general half-off-shell and hence a 
function not only of q and w, but also of the energy and momentum of the off-shell struck nucleon, or equivalently of 
p and £. In the models being considered in the present study the struck nucleon is in fact on-shcU and so TZ^ becomes 
simply the longitudinal response of a moving free nucleon: 

nL = ^{Gl{T) + W2{r)x'} , (2) 

where 

W,ir) = -^[Glir)+rGlir)] , (3) 

expressed in terms of the electric Ge and magnetic Gm Sachs form factors of the nucleon and x — (p/inN) sm9, 
where 9 is the angle between p and q. As in CDM in Eqs. ^ and ^ the dimensionless variables 

A = -^ and T = K^ -X^ (4) 



2mN ' 2mN 

are employed. 

Equation ([T]) connects the semi-inclusive (e, e'N) reaction with the inclusive (e, e') process assuming that the 
outgoing nucleon no longer interacts with the residual {A — 1) nucleus (absence of final-state interactions). That 
equation expresses the assumption that the inclusive cross section is to be obtained by integrating the semi-inclusive 
cross section, summing over struck protons and neutrons. The boundaries of the integration domain E in the {£,p) 
plane are found through the energy conservation relation (see CDM) and read 



£^ {p)=uj-Es+ niN - \l [p ± g)2 + m^ - T^^ " , (5) 

where f + < £ < £^ and £^ is bounded according to 



0<£- <u;-E,- ^q^ + (Ma-i + m^)2 - (Ma-i + m^) . (6) 



In the above 

Es = MA-i+mN~MA (7) 

is the separation energy of a nucleon from the ground state of the nucleus A (of mass Ma) leaving the residual {A— 1) 
nucleus (of mass Ma-i) in its ground state as well and 



4^-'^ = Jp^ + MI_,-Ma-i (8) 



is the very small recoiling nucleus kinetic energy, which is often neglected in Eq. ([5]) . Careful investigations of Eq. ([T]) 
have shown that the singie-nucleon response function is quite gently varying over the integration domain S in the 
(f ,p) plane, particularly if one limits the focus only to regions where the spectral function S{p,£) plays a significant 
role (see CDM); accordingly it can be extracted from the integral in Eq. ^ and divided out, thereby incurring only 
typically few percent errors. This is the point at which the factorization discussed in the Introduction is invoked. 

Finally, in this study we shall confine ourselves to dealing with infinite, homogeneous systems, the simplest among 
them being the RFG model in which the dynamics are controlled by just one parameter, the Fermi momentum kp- 
To attain superscaling it then turns out to be convenient to multiply Eq. ([T]) by a further factor that, in the RFG 
case, is proportional to the momentum transfer and to kp. Indeed, as we shall see, this leads to a function which 
loses any dependence on both kp and q, namely, one has superscaling for the RFG in the non-Pauli-blocked regime. 
It remains to be seen what happens in the BCS model, i.e., in the presence of correlations. 

III. THE SPECTRAL FUNCTION 

As discussed in the previous section the key tool for exploring superscaling is the nuclear spectral function S{p, £), 
defined as follows 

S{p, £) -< -^oWl^SiS -{H- i?(^~^)))apT|*o > -j^ (9) 

in terms of matrix elements of contributions involving the standard fermion creation and annihilation operators. In 
the above |^o > and H are the ground state and hamiltonian of the system, respectively, and Va is the (large) volume 
enclosing the system. 

A. The RFG model 

The RFG spectral function is easily computed when one recalls (see CDM) that it has a negative separation energy 

{Es)rpg = -Tp , (10) 



Tp = y^kp + mj^ — 771 AT being the Fermi kinetic energy. One obtains 

A(27r)3 ' 



5«^«(p,£)= J2 9{kp-p)6{£-Tp + T,)-^, (11) 



spin,isospin 



where the nucleon kinetic energy Tp — ^p^ + ttt,^ — tojv has been introduced, and observes that it is nonzero in the 
(£,p) plane only on the curve 



E = ^kl + ml-^p^+ml (12) 

and is normalized to unity, as is easy to check by integrating Eq. (jlip over p and the non- negative variable £. An 
integration only over £ yields the system's momentum distribution normalized to the Fermi sphere, namely 

df5«^«(p,f) = ^%-3^. (13) 

3^'^P 

Inserting Eq. (jlip into Eq. ^ , after having pulled out and divided by the singie-nucleon response TZl , one obtains 

Fiq,c.) = ?^ // dpd£^S^^^{p,£) , (14) 



where the domain T,iiPG is defined by the curves in Eq. ([5]) using Eq. (fTO]) . Performing the integration, which can be 
done analyticaUy, and dividing by the factor 

1 / rriN \ f 2mNTF 



(note that the quantity 2mi\[TF/kp is very close to unity), one ends up with the REG superscahng function [l| 

fW = 1(1- V^') e (1 - v^^) , (16) 

where 

^=^ , '-^ (17) 

v?^ y'(l + A)t + k^t{1 + t) 

with S^F — TF/rriN. Here ip is the scahng variable specifically obtained within the context of the REG, although it 
has also been widely employed in more general analyses. 

B. The BCS-inspired model 

In the spirit of BCS theory we assume that the ground state is now given by 

\BCS>=l[{uk+Vkai^Ml^^)\0> , (18) 

k 

|0 > being the true vacuum, with 

|wfeP + |wfeP = l, (19) 

to insure the proper normalization of the state, namely 

<BCS\BCS>^l[{\uk\^ + \vk\^) = l. (20) 



With the assumption in Eq. (|18p we have a covariant approximation to the nuclear ground-state wavefunction. We 
have required that the added pairs always occur with back-to-back momenta (hence the net linear momentum of the 
system in its rest frame is still zero) and with opposite helicities (hence the net spin of the ground state is zero; we 
consider only even-even nuclei for simplicity). The creation operators add particles with relativistic on-shell spinors 
u{k, A). 

The state \BCS > does not correspond to a fixed number of particles; in fact it is associated with the spontaneous 
breaking of the C/(l) symmetry in the system. Therefore, it is not an eigenstate of the operator 

n{k) = '^al^aks, (21) 

s 

although its expectation value can be computed, yielding 

n{k) =< BCS\n(k)\BCS >= \vk\^ . (22) 

We can thus require the particle number A to be conserved on the average, which implies the condition 

J2\vk\'^A. (23) 

Concerning the energy, we view our system as being constructed in terms of independent quasi-particles, which 
indeed is the case for the state in Eq. ([T5|) in leading order, writing accordingly 

Ebcs ^< BCS\H\BCS >^< BCS\Y,Ekalaks\BCS >^Y.^k\vk\^ ■ (^4) 

ks ks 



For sake of simplicity in the following we make the assumption 



Ek = \Jml + fc2 , (25) 

which reflects the fact that our present goal is to obtain a reasonable variational wave function (the BCS one), not for 
the purpose of getting the best energy, but rather to provide a way to evaluate the role of correlations on superscaling. 
We now face the question of how to compute the wave functions of the daughter nucleus. For this we assume 

\D{p)>^-^aj>^\BCS'>, (26) 

where Mp is a normalization factor, and 

\BCS'>=\l[u'^{v)+v'^{p)al^al^^]\Q> , (27) 

k 

where the coefficients u'i.{p) and w^(p) are in general different from those in Eq. (fTS]) . but still obey Eq. (fT9|) . This 
point is of crucial relevance for our model, as we shall see below. Note that the daughter nucleus states are labeled 
by the momentum, p. For the normalization factor in Eq. (j26p we get 

Afp^K{p)\' (28) 

and the expectation value of the operator in Eq. (|2ip (the system's momentum distribution) turns out to be (see 
Appendix [X| 

nD^p^{k)=<D{p)\J2alaks\D{p)>=\v',,{p)\\l-S,p) . (29) 

S 

Clearly the condition on the particle number will now read 

E \^'Mf - E i<-(p)i' - K(p)\' = ^ - 1 > (30) 

k^p k 

valid for all of the daughter states, namely for all p. Taking into account Eq. PH]l . this yields the relation 

k k 

to be exploited below. 

The energy of the daughter nucleus in our model is 

EDip) =< D{p)\H\D{p) >= E^fcli'Up)!' , (32) 

k^p 

which, using Eqs. ([24| and ([3T|) . can be recast as follows: 

Ed{p) = {Ebcs - ruN) - Tpl^pl' - E^^tl^'^l' " Wk(P)\''] ■ (33) 

k^p 



We can then proceed to compute the daughter nucleus spectral function 



,(^-1) 



Va 



S'^'^'^p^S) = |< D{p)\ap^\BCS >\'5[S - [Enip) - i^^ '^jj J^ , (34) 

where Eq is the energy E£,rp\ of the daughter nucleus evaluated at that value of p where it reaches its minimum, 

to be referred to as kp in the BCS model. Hence we have 

= Tp\v,,\'~Tp\vp\'+ E n[\v,\'-K{kp)\']~Y.T,[\v,\'~\v'Mn. (35) 

k^kp k^p 



Later we shall show that the last two terms on the right-hand side of the above provide only small corrections and 
so can be dropped. We thus see that our model provides a nice variation of the RFG expression in Eq. (H^); see later 
where results are displayed. The matrix element in Eq. ([34|) can be straightforwardly computed (see Appendix [B| 
and one obtains 



Mp{p) = \< D{p)\ap^\BCS >\ 



Vp\ L ^, W[u'k{p)*Uk+v'k{p)*Vk] 

V ^p(p) feU 



(36) 



Thus we end up with the expression 



^BCS 



ip,£) = \Mp{p)\^S[£-£ip)] 



Va 



A{27ry 



\v,\H[£-£{p)] 



Y[[u'f,{pyuk + v'kipyvk] 



Va 



A{2TrY 



(37) 



for the daughter nucleus spectral function. 

Finally, in order to calculate the longitudinal response function in Eq. ^ what remains to be specified is the 
integration region S, which in turn requires knowledge of the separation energy in Eq. ([7]). In the present model the 
latter turns out to be 






(38) 



The question we are now faced with is how to fix the coefficients Uk, Vk, u'kip) ^^d v'j^{p) entering in the wave 
functions of the initial ground state and final daughter nucleus state. As a first step in this direction let us recall the 
normalization constraints that these coefficients have to fulfill and deduce expressions for them in the thermodynamic 
linnt A — > oo, Va — > oo, A/Va — PA- For the initial nucleus one has 



lim 



-E 



Va 



Vk\ 



IfI"'"|"='"^ 



PA 



(39) 



while for the daughter nucleus one obtains 
1 



lim : 



Va- 



Y.\'^'k(p)?-WM\' 



d^k 2 f Va 

\v {k;p)\ - PA- 



(27r) 



'Va- 



PA-i]\v'{p;p)\' 



lim — = pA-i 

Va-1 



(40) 



Note that if Eq. ([M)) is valid then the spectral function in Eq. ([57)1 is correctly normalized, namely 
/ (Ppd£S'^^^{p, £) = 1, provided that Ylk^pWk^PY^k + <(p)*t'fe] = 1- 

The next constraint for the daughter nucleus wave function's coefficients stems from the stability condition, which 
requires that the energy of the latter fulfill the relation 



dE 



D{p) 



dp 



0, 



(41) 



p=kf 



which, as already mentioned, defines the Fermi momentum kp for the [A — 1) system. Neglecting the last term on 
the right-hand side of Eq. ([33]) this becomes 



which can be recast as follows 

\v{kp)\ 



2 dT{p) 



dp 



I [np)Hp)? 



dv{p) 



= 0, 



P=kf 



+ T{kF)v{kF) , 

p=k^ dp 



+nkF)vikF)'-^ 

p=k, dp 



= 



(42) 



(43) 



p=kf 



Assuming that v is real, this becomes 



v{kf 



y/k' 



l + mjf 



{4^r^„-... 



dv{p) 



dp 



. 



(44) 



p=kf 



In accord with the BCS model and with the physics we intend to explore in the next section, we now choose the 
following three-parameter expression 



v^k) 



.Pik-k) _|_ I 



(45) 



for the V coefficients. Note that when (3 ^ oo one must recover the RFC Substituting Eq. (|45p into Eq. (|44|) we 
obtain 



kp e 



■^ + 1 ^fjJkl+m 



'TV 



kp + mj^ - m^ e 



„0{kp-k) 



(46) 



which can be solved for k in terms of kp and /3, yielding 



k = kp + - log 



I'J''^ 



y fc| + m^ f y fc| + m^ - niN 



-1 



(47) 



The inverse of this equation then gives kp in terms of k and p. 

We use next the normalization condition in Eq. ([39|l to fix the parameter c in Eq. (|45|) . We obtain 



c(/3,fc) = 



n'P'pA 



Li: 



(-e/^fc) 



where 



1 /"^ fc^ 
Li^{z) — — I —r—, r dk 



2 Jo eVz - 1 

is the trilogarithmic function. 

For the daughter nucleus we employ the same form of parametrization as in Eq. (|45p , namely 



(48) 



(49) 



v"{k;p) 



CD 



^Jo[k-~kD{p}] 



(50) 



and require the parameters to fulfill the normalization condition for the A — I system, namely, Eq. (j40p . Assuming 
PA-i — PA = P and Va/Va~i — A/{A — 1) this leads to the equation 



.,««3(-e»^')+^L,3(-e.--«) 



7r2/3' 



P 



CD 



A-1 ePo[p-kD{p)] ^-i 



(51) 



In the thermodynamic limit the right-hand side vanishes and Eq. (j5ip is simply solved by cd = c, Pd = P, kD(j>) = k; 
hence w'^(fc;p) = w^(fc). As anticipated, this implies that the last two terms in Eq. (|35p vanish and that the product 
in Eq. (|37p is unity, so that the spectral function of the daughter nucleus is correctly normalized to unity. It must be 
emphasized that the coefficients v and w' become identical in the thermodynamic limit, but arc different for finite A. 
Hence it is crucial to compute the dispersion relation £(j)) when A is finite and then take the thermodynamic limit. 
This is reminiscent of Koopmans' theorem in Hartree-Fock theory [2J| . 

As far as the parameter /3 is concerned, it clearly controls both the modifications of the momentum distribution 
near the Fermi surface (promotion of pairs due to residual NN interactions, both long- and short-range) and also 
the tail of the momentum distribution due to short-range NN correlations. Indeed, for P very large one recovers the 
familiar ^-distribution, while for smaller and smaller /? more and more particles are pulled out of the Fermi sea and 
produce a significant tail for the momentum distribution at large momenta. The impact of the physics expressed by 
the parameter f3 on the superscaling function is explored in the next section. 




FIG. 1: (Color online) Panel (a): The parameter k, given in Eq. (|47|l . shown as a function of /3. The straight line indicates the 
Fermi momentum kp = 228 MeV/c. Panel (b): The parameter c, given in Eq. (|48|l . plotted versus /3 for pA = kp/{6TT^). In 
both panels the arrow indicates the critical value /3crit = 0.017 c/MeV. 



IV. NUMERICAL RESULTS 



The model is specified by three parameters, k, (3 and c, the last being fixed by the normalization condition in 
Eq. (|48p . The Fermi momentum in this model is determined by the daughter nucleus stability condition in Eq. (|42p . 
Alternatively, as it is somewhat easier to achieve, one can fix the Fermi momentum kp and use Eq. (|47)) to determine 
k. In presenting the results obtained using our model it is then convenient to start by displaying the behaviour of the 
parameters k and c, which are fixed by the physical conditions of normalization and stability, versus (3 for given kp. 
When fc, c and /? are known so are the wave functions of the initial and final nuclei. 

In Fig. 1 the parameters k and c are plotted versus (3. They stay constant (in fact the almost constant value of k 
is quite close to the value kp = 228 MeV/c we use as an input) until a critical value /3crit = 0.017 c/MeV is reached 
where k (c) displays a dramatic decrease (increase). This value corresponds to the change of sign of the logarithmic 
term in Eq. (|47p . namely 



/3c 



2fc, 



Tp{Tp + mN) 



(52) 



For values of f3 lower than /3crit the nuclear momentum distribution becomes very much extended beyond the Fermi 
sphere associated with the input value of kp] we shall return below to discuss this delicate issue. Thus our results 
appear to point to the existence of a narrow domain of [3 around /3crit, below which the system becomes strongly 
disrupted by correlations. This has a strong impact on the structure of the superscaling function, as we shall see 
later. As stated above, for the results shown in the figures we have taken p = kp/{3n^) and kp — 228 MeV/c (a 
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FIG. 2: (Color online) Momentum distribution of the initial state, Eq. (|45)) . evaluated for kp = 228 MeV/c, pA = fc|./(67r^) 
and different values of /3 (in c/MeV) above (a) and below (b) the critical value /3crit- 



value appropriate for ^^C), although our results appear to be quite insensitive to kp as long as it is kept within the 
range appropriate for the physics of atomic nuclei. We display the momentum distribution of the initial nucleus in 
Fig. 2 for a few values of /3 larger (panel a) or smaller (panel h) than ficmt- The progressive development of a tail in 
the momentum distribution is clearly seen in the figure. 

The next important issue to be addressed is to determine where the spectral function is nonzero in the (i?,p) plane. 
The answer is found in Fig. 3 where the support of the spectral functions of the RFG and of our BCS-inspired model 
are displayed and compared for (3—0.1 c/MeV and the same kp used in Fig. 1. Both spectral fmictions of course are 
just (5-functions, but concerning their support two major differences distinguish the two: 

1. In the range of momenta where both exist the excitation spectrum of the daughter system is substantially softer 
than the RFG one; 

2. For missing momenta larger than kp the BCS case, unlike the RFG, continues to display a spectrum (its spectral 
function continues to have a support) which in the thermodynamic limit rises quite suddenly with p until it 
reaches the value £ assumes for vanishing missing momentum, namely 



fmax = Tp\v{kp)Y 



(53) 



This energy, which can be shown to be the upper limit of the daughter spectrum and which cannot be exceeded 
because of the stability condition in Eq. (|4ip . is reached only at p = oo, but over a large span of momenta £ 
remains almost constant, thus corresponding to the situation of an eigenvalue with infinite degeneracy stemming 
from the symmetry C/(l) associated with the particle number conservation. As p is lowered, approaching the 
Fermi surface, the degeneracy is lifted and we face a situation of a spontaneously broken symmetry, reflected in 
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FIG. 3: (Color online) The excitation energy £ computed according to Eq. (|35p neglecting the last two terms for kp ~ 228 
MeV/c, PA ~ k%/(6iY^), P — 0.1 c/MeV. The RFG results are also shown for comparison. 



the structure of our state which contains components of aU possible particle number. This situation is strongly 
reminiscent of superconductivity, where the spontaneous symmetry breaking also occurs in the proximity of the 
Fermi surface (in this case, however, the control parameter is the temperature; in our case it is the momentum 
itself). 

This set of degenerate states has a dramatic impact on the superscaling function of our model, which we obtain using 
the procedure already outlined for the RFG in Section Hill employing the same dividing factor A given in Eq. (jlSp. The 
function / is displayed versus the scaling variable tp in Eq. (|17p for a few values of (3 and q in Fig. 4. For comparison 
the RFG result in Eq. p6|) and the averaged experimental data [3, [2a] are also shown. One sees that to get / for 
large positive -0 we have to integrate in the {£,p) plane in domains (see CDM) encompassing large fractions of those 
degenerate states discussed above. These are thus the cause of the asymmetry of the scaling function with respect to 
"0 = appearing in Fig. 4. For ijj large and negative these states are to a large extent excluded from entering into the 
building up of /. The fact that this effect is more and more pronounced as (3 becomes smaller reflects the impact of 
the tail of the momentum distribution which indeed grows when (3 decreases and, as a consequence, more degenerate 
states participate to build up /. Note that values of (3 around the critical value {j3 ~ 0.01 — 0.02 c/MeV) yield a tail 
which is in qualitative agreement with the experimental data. 

As far as scaling of the first kind is concerned. Fig. 4 shows that this is quickly reached in the vicinity of the QEP, 
although not so to the right and to the left of it. A closer examination of the results (see Figs. 5 and 6, where / is 
plotted on a logarithmic scale for a wider g-range at /3=0.01 c/MeV) shows, however, that also here the BCS model 
does scale, however with an onset reached only for q ~ 1.5 GeV/c, namely for larger momenta than when at the 
QEP where the onset already occurs at about 500 MeV/c. Also from Figs. 5 and 6 it appears that the scaling regime 
is reached faster to the right than to the left of the QEP. Moreover, the asymptotic value for i/; < is approached 
from below, namely the superscaling function grows with q until it reaches its asymptotic value, in contrast with the 
experimental findings. This reflects the fact that our model, although appealingly simple, is not able to account for 
features of this kind. We conjecture that a non-uniform strength function should be associated with the continuum 
spectrum of the daughter nucleus, which in general implies the development of a more elaborate version of our model. 
For instance, a further extension of our model could allow the coefficients v to become complex to account for the 
finite lifetime of the excited states. Note that the same trend of approaching first-kind scaling from below is also found 
in [23, [23, [23| within the framework of the Coherent Density Fluctuation Model where realistic nucleon momentum 
and density distributions are used [2a [23, [23- On the other hand, in relativistic mean-field theory [lO| the approach 
is from above, and thus in better accord with the experimental data. 

Finally, using the present BCS model, we investigate the second-kind scaling behaviour by repeating our calculations 
for different values of the external parameter kp, each one taken to represent a specific nucleus. The results are shown 
in Fig. 7, where it appears that the BCS model breaks scaling of the second kind both at the QEP and for negative -0, 
whereas it fulfills it at high positive values of the scaling variable. However, this is not how the original idea of scaling 
of the second kind was introduced [a, [6| . In that study of experimental data a momentum kA (actually called kp in 
that work, although it was in fact a phenomenological parameter and not necessarily the "true" Fermi momentum 
as it must reflect both initial- and final-state interaction effects) was chosen for each nuclear species and used in the 



11 



=+-, 



exp. I — I — I 

RFG 

/3 =0.1, q =500 

£,=1000 

q =1500 

f3 =0.05, q =500 
'^=1000 
q=l500 

[3=0.02, q=500 

g =1000 

g=1500 

/3=0.01, 9=500 

<? =1000 

q =1500 




FIG. 4: (Color online) The superscaling function / defined in Ea. p4|) plotted versus the scaling variable (|17|) in the RFG model 
and in the present BCS model for three values of q (in MeV/c) and different values of /? (in c/MeV). As usual, kp =228 MeV/c. 
Data are taken from [7ll25ll. 



definition in Eq. (fT7|) of the scaling variable tp and of the dividing factor in Eq. (fT5| , which will now be written as 



falE] [ ^^^nTf 



kA \ q J \ kl 



(54) 



For simplicity, in the present approach the value of kA is chosen in order to have all the corresponding superscaling 
functions coincide at the QEP, thus realizing superscaling at least where the nuclear response is the largest. We 
display the results in Fig. 8, where each curve corresponds to given kp and kA, for q = 1000 MeV/c and /3 = 0.01 
c/MeV. In this case it turns out that scaling of second kind at the QEP can indeed be achieved using the following 
empirical expression 



kA 



1 + a{kQ — kp) 



(55) 



where fco is just a parameter (here ko = kp = 228 MeV/c) chosen to fix the height of the peak and where a = 6 x 10^^ 
c/MeV. That is, kA is found to fall very slightly with increasing kp for realistic values of the Fermi momentum. Over 
much of the range of ip shown in the figure one sees relatively good second-kind scaling, although the results still point 
to a sizable violation of the second kind scaling in the scaling domain (large negative scaling variable). Finally, we 
note that for higher values of P the parameter a is found to decrease, as it should, since the RFG result a = must 
be recovered when /3 — > cxd. The behaviour of a versus the parameter [3 is shown in Fig. 9. 
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FIG. 5: (Color online) The superscaling function / in the negative i/J region plotted for several values of q (in MeV/c) and 
/3 = 0.01 c/MeV. As usual, kp =228 MeV/c. 
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FIG. 6; (Color online) The superscaling function / in the positive t/i region plotted for several values of q (in MeV/c) and 
(3 = 0.01 c/MeV. As usual, kp =228 MeV/c. 



CONCLUSIONS 



In the present study a simple extension of the relativistic Fermi gas model for studies of relatively high-energy 
inclusive electroweak cross sections has been developed. Starting from the RFG in which a degenerate gas of nucleons 
is assumed for the nuclear ground state, in this extension pairs of particles are promoted from below the Fermi surface 
to above, yielding a spectral function and the resulting momentum distribution with Fourier components for all values 
of momentum. In the spirit of the RFG this new model has been constructed in a way that maintains covariance — 
the pairs of particles involved all have spinor products with opposing momenta and hclicitics, thereby producing a 
ground state with net momentum zero and net spin zero. 

Specifically, the model employed is inspired by the BCS approach used in treating superconductivity, both in 
condensed matter physics and in the nuclear physics problem of superconducting nuclei, in which an infinite product 



of operators of the form Up + v^a, 






with quasiparticle weighting factors Up and Vp is assumed. The specific choice 



of these weighting factors characterizes the particular model assumed in the present study, namely, such that the 
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FIG. 7: (Color online) The superscaling function / plotted versus tp for several values of the Fermi momentum kp (in MeV/c) 
at P =0.01 c/MeV and g =1000 MeV/c. 



ground-state momentum distribution has a rounded Fermi surface and a long tail to represent effects from both long- 
and short-range correlations. It is, however, worth recalhng that our state only accounts for correlations between 
pairs of fermions in time-reversed states, namely with opposite momenta, as occur when employing the pairing 
Hamiltonian. Clearly the parametrization used in the present study is only one choice for the weighting factors 
and it is straightforward to employ other more realistic forms when warranted (see for instance [29|, l3d . l3l|)- The 
nuclear particle number and total energy of the ground state are obtained as expectation values using this form and 
constrained to chosen values A and Ma- The final state (the one resulting after the electroweak interaction with the 
ground state has taken place) is chosen to be a plane-wave outgoing nucleon plus a daughter nucleus recoiling in a 
way that conserves momentum. In the rest frame of a general daughter state the particle number and total energy 
may be found and accordingly the daughter ground state may be identified. With these ingredients in hand it is 
then straightforward to obtain the nuclear spectral function and from it to compute the rest of the desired functions, 
namely, the nuclear responses (in the present work only the EM longitudinal response is discussed, although it is 
clearly straightforward to obtain any electroweak responses using the same ideas) , and the superscaling function. 

With respect to the last, the results presented in Section IIVI show that the model is capable of accounting for 
several, but not all, of the features seen in other studies of superscaling, both phenomenologically and within the 
context of alternative modeling. Specifically, scaling of the first kind (independence of q) is observed for momentum 
transfers that are high enough. Near the quasielastic peak the onset of first-kind scaling is quite rapid, whereas away 
from the peak it is somewhat slower, that is, requires a high value of q to occur. As with most models (the exception 
being relativistic mean-field theory and its derivatives) the approach to first-kind scaling is from below, whereas the 
data in the scaling region approach from above. Finally, second-kind scaling (independence of nuclear species) is seen 
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FIG. 8: (Color online) The superscaling function / plotted versus i^} for several values of the Fermi momentum kp (in MeV/c) 
and with A;a devised in such a way that the peaks coincide (see text). 
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to be good although some scaling violations are observed. 

APPENDIX A: MOMENTUM DISTRIBUTION OF THE DAUGHTER NUCLEUS 
The momentum distribution associated with the daughter nucleus reads 
nD{q){k) =< D{q)\ ^ al^aks\D{q) >= 

S 

n K(9)* + v'jMT a-piap^][u'p{q) + v'p{q)al^a^^^y]aksaqi 

p=itq^k 

[u'k{<l) + v'k{q)al^a^_ki][u'q{q) + Wgl'z)^!^-^!]!^ >} = 
P^{l<(9)n«fe('Z)P<0| n K('?)*+<(9)*apT«-Pi]K('^)+<('?)«PT'^-Pi]|0>} 



p=^q^k 



\v',{q)?{^-5kq). (Al) 



APPENDIX B: EVALUATION OF THE MATRIX ELEMENT 



Considering the matrix element 



we observe that 



Mp{q)=<D{q)\ap^\BCS> 
ap^\BCS> = flpT Il("fc + ^fe'*fcT"-fei)l'^ ^ 



n ("fe+«fe4T«-fei 



k=tp 



2pl yup + Vpal^a^_p^^ |0> 



n ("fe + «fe4T«-fei 



fc^p 



^p«-pi|0 > 



-Pi n ("'= + '"'^"iTa-fc J io > 



having defined 



k^p 

Vpal^^\BCS{^p)> 



\BCS{^ p) >= n (uk + Wfe4T«-fci) |0 > ' 

k^p 



the BCS wavefunction without the term p. Similarly 

1 



\D{q) > = 



v'qiqW 



V'q{q) J 

v'qiq)* " 



a,^\BCS' >= ^/-^ al^^lBCS'i^q) > 



and therefore 



Mp{q) = ^pJ^^ < BCS'i^ q)\a^,^al^^\BCS{^ p) > 



(Bl) 



(B2) 



(B3) 



(B4) 
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Spg < BCS'i^ q)\BCSi^p) >-< BCS'i^ q)\al^^a-,^\BCSi^p) > . (B5) 



Now 



Sp, < BCS'i^ q)\BCSi^ p) >= Sp, llKiqyuk + v'kiqYvk] , (B6) 

whereas the second term in Eq. (|B5[) vanishes, since 

a_5x|SC5(^p)>=0 iip^q (B7) 

and 

a_g^\BCS{y^p)>^^Vgal^\BCS{^p,^q)> ii p ^ q, (B8) 

namely 

a_,x \BCSi^ p) >= (,5p, - l)vqal^\BCS{^ P,^q)> (B9) 

and, similarly 

a^p^\BCS'{^q)>={5p,-iyj,{q)al^-\BCS'{^p,^q)> . (BIO) 

Then 

< BCS'i^ q)\a-gial^^\BCS{^p) > 

= {dp, ~ l)<(g)*t;, < SC5'(^ p, ^ g)|apT4^|SC5(^ p, ^ g) >= . (BU) 

Eqs.([B6| and (|BTT|) . when inserted into Eq. (|B5l), finally yield 



/ v' (q) 
Mp{q) = 5p,v,J^^Yl[u',{qrUk + v',{qrVk] . (B12) 
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